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We extend a recent formulation of quantum continuum mechanics [J. Tao et. al, Phys. Rev. 
Lett. 103, 086401 (2009)] to many-body systems subjected to a magnetic field. To accomplish 
this, we propose a modified Lagrangian approach, in which motion of infinitesimal volume elements 
of the system is referred to the "quantum convective motion" that the magnetic field produces 
already in the ground-state of the system. In the linear approximation, this approach results in 
a redefinition of the elastic displacement field u, such that the particle current j contains both an 
electric displacement and a magnetization contribution: j = jo+no<9tU+ V x (jo x u), where no and jo 
are the particle density and the current density of the ground-state and dt is the partial derivative 
with respect to time. In terms of this displacement, we formulate an "elastic approximation" 
analogous to the one proposed in the absence of magnetic field. The resulting equation of motion 
for u is expressed in terms of ground-state properties - the one-particle density matrix and the 
two-particle pair correlation function - and in this form it neatly generalizes the equation obtained 
for vanishing magnetic field. 



I. INTRODUCTION 

Consideration of an hydrodynamical formulation of the 
electron dynamics goes back to the early days of quantum 
mechanics^. In the modern language of time-dependent 
(current-)density- functional theory^ (TD(C)DFT), this 
naturally leads to appealing time-dependent orbital-free 
methods^. In this spirit, in two recent papers^, 
the problem of calculating the linear response of a 
generic quantum many-body systems to an external time- 
dependent potential has been reformulated in the lan- 
guage of quantum continuum mechanics (QCM). In this 
approach, the non-equilibrium state of the system is de- 
scribed in terms of an elastic "displacement field" u(r, t) 
(a function of cartesian coordinates r and time i), such 
that an infinitesimal volume clement of the system that is 
located at point r in the equilibrium state will be located 
at r + u(r, t) in the non-equilibrium state. The particle 
current density is connected to the displacement by the 
relation 



j(r,t) = n (r)d t u(r,t) , 



(1) 



where dt represents a partial derivative with respect to 
time and no(r) is the ground-state density. An "elastic 
approximation" was then introduced, based on the idea 
that the time evolution of the wave function can be de- 
scribed as a geometric deformation of the ground-state 
wave function, the deformation being defined by the dis- 
placement field u. More precisely, the wave function of 
the deformed state IV 1 ! 11 ]) was expressed in terms of the 
ground-state wave function \ipo) in the following manner: 



\ip[u]} = exp 



dv j(r) -u(r) 



(2) 



where j (r) is the canonical current density operator act- 
ing as the generator of differential translations - a dif- 
ferent translation at each point in space. (Here and in 
the following, we set the mass of the particles m = 1). 
Starting from Eq. ([2]), a closed equation of motion for 
u could be derived, which is demonstrably exact for (i) 
one-particle systems, and (ii) many-particle systems sub- 
jected to high-frequency fields. This equation has the 
form 



n (r)9 t 2 u(r,f) 



-n (r)u • V(VPq) 



ST 2 [u] SW 2 [u] 



5u(r, t) <5u(r, i) 
-no(r)Wi(r,t), (3) 



where Vb(r) is a static potential that defines the na- 
ture of the many-body system, V\ (r, t) is the small time- 
dependent potential to which the many-body system re- 
sponds, T 2 [u] is the kinetic energy density of the de- 
formed state @ expanded to second-order in u, and 
W2 [u] is the electron-electron interaction energy density 
of the deformed state, also expanded to second-order in 
u. The functional derivatives that appear on the right 
hand side of Eq. Q are actually linear integro-differential 
operators acting on u. The explicit form of these op- 
erators was derived in Refs. [fjjj?]. The final expressions 
involve only the following ground-state properties: the 
one-particle density matrix and the two-particle correla- 
tion function, which arc quantities that may be computed 
by means of quantum Monte Carlo methods^. 

Quantum continuum mechanics holds great promise as 
a tool for simplifying and streamlining the calculation 
of excitation spectra, particularly in situations in which 
the ground-state correlations are well understood and the 
spectrum is dominated by collective excitations. Another 
interesting possibility is to use quantum continuum me- 
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chanics as a tool for efficiently approximating the density- 
density response function of the non-interacting Kohn- 
Sham system. This possibility has been recently pur- 
sued by Gould and Dobson^ The Kohn-Sham response 
function is then used by these authors, in combination 
with the random phase approximation, to generate more 
accurate exchange-correlation energy functionals, which 
capture dispersion forces between metals and insulators. 
Moreover, the use of quantum continuum mechanics has 
allowed to better understand and simplify the deriva- 
tion of the expression of the high-frequency limit of the 
exact exchange-correlation kernel of TD(C)DFT— . We 
should also mention that, another interesting approach 
that holds promise to overcome present limitations of the 
available approximate KS methods is the semiclassical 
density-matrix time dependent propagation of Rcfsi 11 ' 12 . 
The theory of Refs. BU 

was based on the assumption 
that the ground-state of the many-body system is time- 
reversal invariant and has no spin-orbit coupling, so that 
the ground-state current-density jo(r) = 0. This paper 
is concerned with the extension of that theory to many- 
electron systems in the presence of a static magnetic field 
Bo(r) = V x Ao(r), where Ao(r) is a static vector poten- 
tial. Thus, we consider a system of N identical particles 
described by the time-dependent Hamiltonian 



N 



1 

2 



HV,+A (r,)) 2 



+ V (r j ) + V 1 (r j ,t) 



(4) 



where W(\r — r'|) is the electron-electron interaction po- 
tential (we set e = h = c = 1). The time-dependent 
many-body wave function \&(ri, . . . , r^r, t) is the solution 
of the Schrodinger equation 

idtV(r 1 ,...,T N ,t) = H^(r 1 ,...,r N ,t) (5) 

with initial condition 



*(n,...,rjv,0) = *o(ri,...,r7v). 



(0) 



There are compelling reasons for working out the gen- 
eralization of QCM to systems subjected to magnetic 
fields. As discussed in Refs. 00, the elastic approxi- 
mation is, in essence, a collective approximation for in- 
homogeneous systems. It condenses the excitation spec- 
trum of the many-electron system into a simpler spec- 
trum of collective excitations, which still carry the exact 
full spectral strength (this is a consequence of the exact- 
ness of the theory in the high-frequency limit). Naturally, 
such an approximation becomes more trustworthy when 
the excitations under study are truly collective, as op- 
posed to incoherent single-particle excitations. But, it 
is well known that a strong magnetic field, by quench- 
ing the kinetic energy of an electronic system, suppresses 
single-particle behavior and promotes collective behav- 
ior. Indeed, one of the first successful theories of the ex- 
citation spectrum of the homogeneous two-dimensional 



electron gas (2DEG) at high magnetic field was based 
on a single-mode approximation very similar to our elas- 
tic approximation^. Nanopattcrncd electronic systems 
at high magnetic field (e.g., quantum dots) also exhibit 
strongly collective behavior. One can, for example, fab- 
ricate a lattice of quantum dots (nanopillars) by chemi- 
cal etching on a 2DEG^ and observe the collective ex- 
citations of the resulting electronic system by Raman 
scattering^. By doing this experiment, new collective 
modes have been recently discovered^, which have no 
counterpart in the homogeneous 2DEG. We believe that 
our continuum mechanics will be useful precisely for a 
microscopic study of these collective modes. 

The generalization of Eq. ([3]) to systems described by 
the Hamiltonian (|4]) is not as straightforward as one 
might initially think. Of course, the presence of the vec- 
tor potential modifies the form of the kinetic energy and 
introduces a Lorentz force term, but this is not the main 
difficulty. The main difficulty arises from the fact that 
the ground-state current no longer vanishes: it has a fi- 
nite expectation value, jo(r). This introduces an ambi- 
guity in the definition of the displacement field. 

The first possibility we explored is to define u in close 
analogy to Eq. (fTJ) such that 



j(r, t) = jo(r) + n (r)<9 t u(r, t) 



(7) 



With this definition, it is possible to derive within stan- 
dard linear response theory, a closed equation of motion 
for u. This will be done in Section II. While this equation 
is formally elegant and gives some insight into the gen- 
eral properties of the solutions, it is very difficult to put 
it in an explicit and therefore useful form. The reason is 
that one needs to calculate the ground-state expectation 
value of complicated commutators: the amount of alge- 
bra involved is formidable. By contrast, in the treatment 
of the zero field case we could rely on a direct calculation 
of the energy of the deformed state © - a comparatively 
simpler task that did not require the evaluation of com- 
plicated commutators. 

Eq. §7$ assumes that the excess current j — jo is en- 
tirely due to the time derivative of the displacement field, 
but, in the presence of a magnetic field, even a time- 
independent displacement can produce an excess current 
(see below). After all, the ground-state is perfectly sta- 
tionary, and yet it does carry a current jo- However, at 
variance with a time-dependent current, the current as- 
sociated with a static deformation must necessarily have 
vanishing divergence in order to satisfy the continuity 
equation. Thus, for example, one must have V • jo = 0. 
Taking into account this condition, we are free to add to 
the right hand side of Eq. ((7|) the curl of a "magnetiza- 
tion field" . In other words, we expect the most general 
form of the relation between current and displacement to 
have the form 

j(r, t) = j (r) + n (r)d t u(r, t) + V x M(r, t) , (8) 

where Ai(r, t) is a functional of u. The divergence of the 
last term on the right hand side is guaranteed to be zero. 
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This is completely analogous to the material current in 
electrodynamics, which is customarily written as the sum 
of the time derivative of the electric polarization and the 
curl of the magnetization^ 

To determine the form of A4(r,t) up to the linear or- 
der in u, we assume that (i) M.(r,t) is a local func- 
tional of the displacement, i.e., it depends on u(r, t) at 
the same point of space and time, (ii) A uniform displace- 
ment u(r) = u must cause the ground-state current to 
be rigidly displaced by jo(r) — > jo(r — u), i.e. we must 
have V x J\A(r,t) = — (u • V)jo(r) for uniform u. This 
condition implies that Al(r) = u x jo(r) and fixes the 
relation between current and displacement in the form 

j(r,f) =jo(r) + n (r)<9 t u(r,t) + Vx [u(r,i) x j (r)] . (9) 

Remarkably, this relation emerges naturally from the 
Lagrangian formulation of the problem, which we present 
in Section III. We remind the reader that in Refs.dQ the 
elastic approximation was derived via a transformation 
to a local non-inertial reference frame - the so-called co- 
moving frame - in which the density is constant and the 
current density is zero. To achieve these conditions, the 
displacement field u, which defines the transformation 
to the co-moving frame, was related to the current via 
Eq. (TTJ. In the present situation, the physically "nat- 
ural" requirement for the co-moving frame is that the 
density remain constant and the current density remain 
equal to the ground-state value jo- In other words, an 
observer "riding" on a volume element should not detect 
any change in the density or the current density. We 
found that this requirement determines the relation be- 
tween the displacement field and the current density in 
the form of Eq. (|9]) just as we found from the heuristic 
argument given above. 

Adopting the mentioned special co-moving frame is the 
crucial insight that allows us to arrive at an explicit equa- 
tion of motion for u in the presence of a magnetic field. 
Throughout the paper, we will consider a general situ- 
ation of non-collinear and non-uniform magnetic fields 
(neglecting spin-degrees of freedom). The equation of 
motion in the practically most common case of a uniform 
magnetic field is obtained from the equation of motion 
§3§ by the following simple replacements: 

1. Replace the time derivative dt by the "convective 
derivative" 

A = dt + v • V (10) 
where v = jo/no- 

2. Include on the right hand side a "Lorentz force 
term" 

AuxB . (11) 

3. Calculate T2[u] and respectively from the ex- 
pectation values of the kinetic energy and electron- 
electron interaction energy operators evaluated in 



the deformed ground-state wave function defined, 
just as in Eq. @. The electron-electron interac- 
tion energy term remains formally unaffected, while 
the kinetic energy term is modified by the inclu- 
sion of an effective vector potential accounting for 
the external static vector potential and the cor- 
responding contribution of the convective motion 
in the ground-state (remarkably, this modification 
vanishes for one-electron systems). 

The paper is organized as follows: 

In Section [TT] we derive an approximate equation of 
motion for the conventional displacement field from the 
high-frequency expansion of the exact current-current re- 
sponse function. We show that this approximation is ex- 
act for one-particle systems and discuss the difficulties 
arising when one attempts to explicitly evaluate the for- 
mal expressions appearing in this equation. 

In Section Mil we address the difficulties discussed in 
Section|TT]by resorting to a non-standard Lagrangian for- 
mulation appropriate for systems in magnetic field. To 
this end, we formulate the quantum many-body dynam- 
ics in a special co-moving frame, such that both the den- 
sity and the current density retain their initial (ground- 
state) values at all times. We arrive at an exact, but still 
not explicit equation of motion for the displacement field 
in terms of the Hamiltonian in the stress tensor of the 
co-moving frame. 

In Section IIVI we introduce the elastic approximation 
and obtain a closed, fully nonlinear equation of motion 
for the displacement field in this approximation. 

Finally, in Section [V] we linearize the elastic equa- 
tion of motion and obtain an explicit, linear equation 
of motion for the displacement field. It is shown that 
the displacement field obtained in this manner is not the 
conventional one, but is related to the current density by 
Eq. ©. 

Section I VII concludes the paper with a summary of the 
main results. 



II. DERIVATION FROM LINEAR RESPONSE 
THEORY 

A formally exact equation of motion for the cur- 
rent density response of a many-particle system in 
the linear response regime can be easily derived from 
standard linear response theory 19 ! 20 To this end, it 
is convenient to replace the external potential V\ (r, t) 
by a time-dependent vector potential Ai(r, t) — 
— J_ VV(r, t')dt' , which is physically equivalent since 
it gives rise to the same time-dependent electric field and 
no time-dependent magnetic field. Assuming further that 
the time-dependence of the external field is periodic with 
angular frequency oj and switched on adiabatically at 
t = — oo with the system initially in the ground state 
\ipo), we obtain the standard result for the Fourier com- 
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ponent of the current density at frequency u>: 

ji, ll (r,u) = f dr' Xfw (r,r',u)^^l, (12) 
J iu 

where X jUt /(r, r', w) is the current-current response func- 
tion (Einstein summation convention is used throughout 
the paper) 

X^(t,t',u) = n (r)5(r-r') 



and 



9„Vi(r') 



* / dte^(Vo|[j(r,t)J(r',0)]|Vo) 
/o 

(13) 

is the Fourier component of Ai(r,i) at fre- 



quency ui. In the above equation (flU| j(r,i) is the 
current-density operator defined as 

1 N 

fr) = 2 E {H^i + Ao(*i)], 5(r - f,)} , (14) 

i=i 

where {A, B} = y4B + £?A is the anticommutator. This 
operator evolves in time under the unperturbed Hamil- 
tonian, i.e., the Hamiltonian Q without the V± term. 
The first term on the right hand side of Eq. (|12[) . is 
the so-called diamagnetic response, which involves only 
the ground-state density n (r). A formal inversion of 
Eq. (|12j) yields immediately an equation of motion for 
the current density: 

iuj / dr'[x _1 ] MI ,(r,r',w)ji iiy (r',w) = d^V^r,^) . (15) 



Obviously, this result is purely formal, since we have 
no way to exactly calculate the current-current response 
function of a many-body system, let alone invert it. How- 
ever Eq. (|15[) can serve as a convenient starting point for 
constructing approximate theories. Following the ideas 
proposed in Refs£i£ we consider the high frequency limit 
of Eq. (|15p and then interpret it as an approximate equa- 
tion of motion for the induced current density. 

We start by observing that, at high frequency, the 
current-current response function has the well-known ex- 



pansior 



19,20 



Xvu,(r,r',u}) = n (r)5(r - r')S^ - i 



*Wr,r') 



where 



and 



M^(r,r') = -(M[HoMr)]Mr')]\ih) 



(16) 



(17) 



(18) 



A few words of caution should be added at this point . 
The expansion written for x m Eq. (|16[) is done under the 



assumption that the real-space current-current response 
function has a regular Taylor expansion in inverse pow- 
ers of at high-frequency. However, this is not always 
the case. If the initial (ground) state is not sufficiently 
smooth, the unboundedness of the kinetic energy opera- 
tor may cause Xm"( w ) ^° develop a non-analytic behav- 
ior (e. g. fractional powers of ui) at high frequency^. 
Therefore, strictly speaking, the expansion of Eq. (|16[) 
assumes that proper smoothness conditions arc imposed 
on the initial state of the system. Fortunately, after in- 
verting Eq. (|16j) to get the high frequency expansion of 
the operator X ~l entering the left hand side of Eq. (fT5|) , 
all smoothness restriction can be relaxed. It turns out 
that such obtained high frequency form of the inverse 
response function is generally valid. We prove this in 
Sees. IV and V by redcriving the equation of motion for 
the current density via the Lagrangian frame formalism, 
which can be viewed as a direct construction of a high 
frequency / short time limit of Xal ■ In Appendix B we ex- 
plicitly demonstrate that for the one-particle system the 
above formal inversion procedure indeed yields the exact 
form of the inverse current response function. A possible 
nonanalytic behavior in ui of the response is correctly re- 
covered in our theory because at the level of Eq. (|15l) it 



is encoded in the space part of the operator xZv acting 
on the current density. 

Thus, we invert Eq. (|16p and plug the result into 
Eq. ((151). This yields 



u + iuj J dr ' B t^( v > v 

J dr'/C^(r,r') 



no(r') 
no(r') 



d»Vi (19) 



where 



/C M „(r,r') = M^(r,r') 

+ J rfr% 7 (r,r'')^pyB 7 ,(r'',rO.(20) 

Evaluation of the current-current commutator in 
Eq. (fT7|) is relatively straightforward and yields 



*Wr,rO =io ) ,(r0^f^ + io 1 ,(r)-^ (r '- r) 



dr v 

+ e ftI/7 S , 7 (r)n (r)(J(r - r') 



dr u 



(21) 



where jo(r) is the ground-state current. Notice that 
E^„(r, r') is imaginary and antisymmetric. 

The evaluation of the double-commutator K.^{v : r'), 
while in principle equally straightforward, is in practice 
an extremely cumbersome task, leading to very compli- 
cated expression to which we are not able to attach any 
transparent physical meaning. We do not undertake this 
task here, since in the next section we will develop an al- 
ternative approach, which leads more directly to a phys- 
ically meaningful equation of motion. 

Even without knowing the explicit form of /C we can 
learn something from the form of Eq. (|19|) . To begin 
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with, we may express the response of the current in terms 
of the displacement defined in the conventional way, i.e. 
according to Eq. 0. By doing so, we obtain 

J 1 n u + iu)Bu — Ku = n Wi , (22) 

where we have adopted a compact notation in which 

O u = Jdr' CV(r,r'K(r') . (23) 

It is easy to prove that M must be a positive-definite 
operator, since, by definition, it is related to the second- 
order term in the expansion of the energy of the distorted 
ground-state 10 in powers of u. 

Setting the right-hand side to zero in Eq. (|22|) . we rec- 
ognize a generalized eigenvalue problem of the form (for 
more details see Appendix A) 

cj 2 n u + icjS-u-7C-u = 0. (24) 

Under the stronger assumption that not only M, but 
also 7C is positive definite, one can show (see Appendix 
A) that the generalized eigenvalue problem in Eq. (|24|) 
has real solutions ui whose modes satisfy a generalized 
orthogonality relation 

(u B ,iB-u A ) + (u> A + oj b ) (u B ,n Q u A ) = , u> A ^ u> B 

(25) 

with the scalar product defined as follows 

(u b ,ua) = J dr UB iM (r)uA, M (r) . (26) 

The eigenvalues are naturally interpreted as (approx- 
imate) excitation energies and the corresponding eigen- 
vectors are matrix elements of the current density oper- 
ator between the ground-state and the excited state in 
question (see Ref. and Appendix B in the present pa- 
per). Thus, we see that the assumption of positivity of 
JC is, in practice, equivalent to the expectation that our 
approximation does not lead to spurious instabilities (i.e. 
complex excitation energies). 

A natural question at this point is: how reliable our 
approximate equation of motion (|22l) will be in situations 
other than the very high frequency limit? The answer to 
this question has already been discussed extensively in 
Ref. l7|, so we only summarize the main points adding 
the required modifications where needed. 

First of all, the fact that the proposed current-current 
response function has the exact high-frequency behavior 
up to order 1/uj 2 implies that the first moment of the 
spectral function dui ui^mx^u (r, r', d) ~ M^r, r') is 
exactly reproduced. Thus, while the energy of individual 
excitations may be misrepresented, the first moment of 
the current-density fluctuation excitation spectrum (or, 
equivalently, the third moment of the density fluctuation 
excitation spectrum) is correctly reproduced. 

Second, our equation assumes an exact and explicit 
form for one-particle systems. This is rigorously proved 



in Appendix B and Appendix C. In particular, in Ap- 
pendix B we show that % _1 (w) ji = A+iB/ui + C/ui 2 for 
any frequency u>, where A, B, and C are coefficients ex- 
pressible in terms of the ground-state density and current 
density. The coefficient A is related to the diamagnetic 
response, while B and C may be related to the coeffi- 
cients of lu^ 1 and u)~ 2 in the high-frequency expansion 
of X- 

The fact that the form derived by the above proce- 
dure becomes exact in one-particle systems has a deeper 
physical significance. Recall that in Rcfs. 00 the equa- 
tion for the displacement was derived by describing the 
dynamics in a co-moving reference frame, defined as an 
non-inertial frame in which the density remains constant 
and the current density is always zero. So the high- 
frequency approximation could be restated in terms of 
an "anti-adiabatic approximation" or "elastic approxi- 
mation" introduced as the assumption that not only the 
density and the current, but the wave function itself re- 
mains unchanged in the co-moving reference frame. This 
assumption is correct for one-particle systems, because in 
such systems the density and the current density uniquely 
determine the wave function, up to a trivial phase factor. 

We will show that this approach - transformation to 
a co-moving frame and the requirement of a stationary 
wave function in that frame - also works in the presence 
of a magnetic field, but with an important difference: the 
current density in the co-moving frame will not be zero, 
but will be equal to the current density in the ground- 
state. Because of this redefinition the relation between 
current response and displacement field will be drasti- 
cally changed: Eq. (fTJ) will be replaced by Eq. ([8]). How- 
ever, as a result of this modification, we will be able to 
obtain an explicit equation of motion for the new u, by- 
passing the need to calculate cumbersome commutators. 
The final form of the equation of motion will be a natural 
extension of the one derived in Rcfs. in the absence 
of a magnetic field. 



III. LAGRANGIAN FORMULATION 

The key physical idea of the Lagrangian formulation 
in application to many-body dynamics is to separate the 
convectivc motion of the electron fluid from the motion of 
electrons relative to the convectivc flow. The former type 
of motion is characterized by the trajectories of infinites- 
imal volume elements, while the latter is described by a 
many-body Schrodinger equation in a local non-inertial 
reference frame attached to those elements 2 ^. 

In the standard Lagrangian formalism,— the convec- 
tive dynamics of the system is described by a set of 
trajectories r(£,t), which represent the motion of an in- 
finitesimal volume elements initially located at These 
trajectories satisfy the first-order differential equation 
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with initial condition r(£, 0) = Notice that we are not 
even making the linear response approximation at this 
stage: the quantity n(r(£, t), t) is the actual particle den- 
sity along the trajectory and reduces to the ground-state 
density Jio(£) only at the initial time. The displacement 
held u(£, t) is dehned via the relation r(£, t) = £ + u(£, t). 

This standard definition is, however, inconvenient if 
the initial stationary/ground state already carries a non- 
zero current-density j'q = uqVq , which is the case in 
a magnetized electronic system. Indeed, even in the 
absence of the external driving field the ground state 
flow drags the volume elements and produces a time- 
dependent displacement in a system that, physically, 
should be considered as stationary and undeformed. It is 
natural to try and exclude the equilibrium convective mo- 
tion from the trajectory function, i.e., to consider a vol- 
ume element as moving only when it moves relatively to 
the ground-state flow. To accomplish this in the presence 
of the magnetic field we adopt a modified Lagrangian de- 
scription, where the motion of the volume elements is de- 
scribed relatively to the ground state flow. Thus, when 
the system is in the ground-state the displacement van- 
ishes and the volume elements are regarded as stationary. 
An observer "riding on" these volume elements does not 
move at all (relatively to the laboratory) and sees the 
current density jo- Now, when the system is out of equi- 
librium a finite displacement appears and the volume el- 
ements begin to move. However, to an observer attached 
to the material elements the current density still appears 
to be jo- Just as in the standard Lagrangian description 
an observer riding on the volume element sees stationary 
density Uq and zero current density, in the present de- 
scription an observer riding on the volume element sees 
stationary density no and current density jo- 

The task at hand is now to implement the transforma- 
tion that makes the density and the current density equal 
to their ground-state values at all times. This is done in 
two steps. First we learn how to describe the quantum 
many-body dynamics in a generic non-inertial reference 
frame. Then we fix the reference frame from the condi- 
tion that density and current density retain their ground- 
state values at all times. We carry out the above steps 
without assuming that the displacement is small, i.e., in 
a fully non linear way. The linearization of the equations 
of motion will be carried out only at the very end of the 
next Section. 



A. Transformation to a non-inertial reference 
frame 

Let us consider a local reference frame moving along a 
prescribed trajectory r(£,t). The form of the quantum 
many-body dynamics in such a frame is worked out in 
full detail in Refs. [25I - 42TI . Here we present only the key 
results that are needed for our purposes. 

First of all, the transformed many-body wave function 



has the form 



JV 



where 



x ^>{v{^t),...,v{i N ,t),t), (28) 
dr a (£,t) dr a {£,t) 



(29) 



is the metric tensor induced by the (non-singular) trans- 
formation r — > £, y/g = \J det g^ is the Jacobian of the 
same transformation and 



dt 



1 



-(r(£,i)) 2 -r(£,i)A(r(^),i) 



(30) 



is the classical action of a particle moving along the as- 
signed trajectory. The factor Ylf =1 g^(£j,t) in Eq. (|28|) 
preserves the standard normalization of the wave func- 
tion (\E , |\E r ) = 1 after a non- volume-preserving transfor- 
mation of coordinates. Equation (|28[) is a generaliza- 
tion of the transformation to a homogeneously acceler- 
ated frame, which is used, for example, in the proofs of 
a harmonic potential theorem^—. 

The transformed wave function "f^i, . . . , t) satis- 
fies the Schrodinger equation 

id t $(Z 1 ,...,Z N ,t)=H\g llv ,AMZ 1 ,...,£ N ,t) (31) 
with the transformed Hamiltonian 



N 



H[g^A]=Y,9j 1 K : 



9 3 9, 



3 = 1 



(32) 



where 



dr v dr v 
gf^oAr(U))- W r^r(ttU) 

dpS*(t,t), (33) 



and l£ k £j is the distance between jth and fcth particles 
in the non-inertial frame (i.e., the length of geodesic con- 
necting points £j and £k in the space with metric g^ v ). 
The metric tensor g^ }V and "effective" vector potential 
A in Eq. (|33|) describe the combined action of external 
forces and inertial forces in the local non-inertial frame. 

The transformed densities are obtained from the re- 
duced one-particle density matrix 

N 

through the formulas: 

=7(C, €,*), (35) 
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^')7(€,€',t) 



(37) 



where 

i 

jpA&t) = 9 J} m tu ( d t» 

is the (covariant) paramagnetic current. It is extremely 
important to carefully keep track of covariant (lower in- 
dices) and contravariant (upper indices) components of 
vectors and tensors. The two types of components are 
connected by the metric tensor via the standard formu- 
las 



(38) 



where g^ v is the inverse of g^ v . Armed with these defi- 
nitions, one can readily verify that the expressions of the 
densities in the non-inertial frame are related to the ones 
in the laboratory frame by the relations 

»»(€,*) = Vgn(r(£,t),t), (39) 

~ne,t) = v^^[i M (r(€,*),t) 

- n(r(t,t),t)r»(T($,t),t)] . (40) 

The local conservation laws for the transformed many- 
particle problem read as follows 



d t n + d^f = 



(41) 



and 



d t j^f{d^A^-d^A v )-nd t A^+^D v P; = 0, (42) 

where P^ is a mixed component of the stress tensor, 
which, in Ref. [25| was proved to have the form 



P^(£,i) = -— 



5H[g af} ,A] 



(43) 



where the functional derivative with respect to g^ is to 
be taken fixed A~ 

The quantity D U P^ in Eq. (|43|) stands for the covariant 
divergence of the stress tensor and is explicitly given by 



A rather lengthy calculation shows that 



dr v ~ 




(44) 



(45) 



where the functional derivative with respect to r v is to 
be taken at constant A. It is important to notice that 
in the above equation, the transformed hamiltonian H 
- a functional of g^ v and A [see Eq. (|32|) ] - is actually 
treated as a functional of the trajectory r M and A. This 
is permissible because g^ itself is a functional of r M - see 
Eq. (USD. 

We will make heavy use of this identity in the next 
section. Now let us learn how to fix our reference frame 
so that the density and current density become constants 
of the motion. 



B. Fixing the reference frame 

As discussed in the introduction to this section, our 
goal is to separate the convective motion of the electron 
fluid from the motion of electrons relative to the con- 
vective flow. This is achieved by following the volume 
elements along their trajectories, so that the density and 
the current density becomes stationary, with values equal 
to the initial one. Mathematically, this translates to the 
single condition 



(46) 



It is important to notice that our condition must be 
imposed on the contravariant components of the cur- 
rent density, rather than on the covariant components j M . 
The two choices are, at first sight, equally plausible, but 
not equivalent, since the metric tensor c/ M „, which con- 
nects covariant and contravariant components of vectors, 
is time-dependent. What forces our choice is the fact 
that the contravariant component of the current density 
appears naturally in the continuity equation (|41[) . Thus, 
with this choice, the continuity equation (|41j) guarantees 
that we get 



»(€.*) =n(€.* = 0) = no(€) 



(47) 



as soon as Eq. (|46j) is satisfied. 

Given Eq. (|4"6"j) and Eq. (14"T1) . the effective potential 
acting on the system is determined through Eq. (|36[) as 
follows 



Mi) ' 



(48) 



where Vq(£) = j£ (£)/ n o(l) • Eq. (gSJ) is a gauge-fixing 
condition which determines the effective vector potential 
A. and, thus, a particular non-inertial frame. 



C. Equation of motion for the Lagrangian 
trajectory 

After fixing the reference frame the system is com- 
pletely characterized by two dynamical variables - the 
trajectory of volume elements r^(£,t), and the trans- 
formed many-body wave function ^(t) that describes mi- 
croscopic motion relative to the convective flow. The 
wave function ^(t) satisfies the Schrodinger equation 
pil) . In order to obtain a physicaly trasparent form of the 
equation of motion for the Lagrangian trajectory r^(£, t) 
we differentiate Eq. (f33|) with respect to time and we get 

f + ^Fo(r) + [r x B (r)]„ + ^d t A v = -0„Vi(r,t) . 

(49) 

In this equation we insert the value of d t A v determined 
through the local momentum balance equation [Eq. (I42[) ]. 
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IV. ELASTIC APPROXIMATION 



nodtAn = d t (flvio) - Jo ~ 9^A V ) 

+ yfgD v P; , (50) 



where we have made use of the conditions (|46|) and (|47|) 
and lowered the index of the contravariant current j„ = 
9nuj u through the action of the metric tensor g^ u . And 
now, on the right hand side of Eq. (|50|) we plug in the ex- 
pression ([33| for Afj, in the combination {d^A^ — d^Av) 
(notice that the contribution of the classical action term 
vanishes, since the curl of a gradient is zero), and make 
use of the identity (|45|) for ^/gD^P^. The result of these 
manipulations is 



2vZd, 



1 

n 



where 









(' 








\Sr v J 








A. 




D t = d 


t + 



(D t r) x Bo^ + ^oW 

*)= -dpVi(T,t) , (51) 



(52) 



is a "convective time derivative" , which takes into ac- 
count the ground-state flow. 

In view of the last result, it is desirable to change all 
the time derivatives in Eq. (|51[) to convective ones. To 
this end, we note that 



r» + 1vld v r» = (5 t 2 + 2v v d v dt) r» 
and "complete the square" to get 

1 5Wo{g^} 
no (5r M 



(53) 



(54) 



where 

W { 9flu ] 



d€no(Off^tt.*K«K(€)- (55) 



Using Eq. ([54]) in Eq. (|5T1) . we obtain 
D?r" + [(Ar)xB o ] M + M Vo(r) 





<5Wn 



= -^Vi(r,t) 
(56) 



We stress that the latter equation, Eq. is an exaci 
nonlinear equation for the trajectory of current-carrying 
material elements in the quantum many-body system, 
provided ^(t) entering the stress force is the solution to 
the Schrddinger equation (|31[) . 



We now introduce our elastic approximation by setting 

= *o(€), (57) 

It is clear that this approximation is consistent with the 
requirement of stationary transformed densities, but can- 
not be exact in general, since the correct form of ^ is 
determined by Eq. (f3"T|). 

The approximation in Eq. (|57p applies to short-time 
dynamics, or to fast- driving fields such that there is no 
time for the wave function to adjust to minimize the en- 
ergy in the presence of fast- varying external conditions. 
Here, we may recognize the high-frequency approximation 
discussed in Section [XXJ The same approximation may be 
characterized as an anti-adiabatic approximation in the 
following sense. The wave function in the Eulerian (iner- 
tial) frame is obtained as an instantaneous deformation 
of the initial wave function. While the initial wave func- 
tion minimizes the energy (i.e., it is the ground-state), it 
is clear that the deformed wave function does not mini- 
mize, at any given instant, the energy of the system. The 
deformation is similar to the change in the shape of an 
clastic body. Hence, the term elastic approximation. 

Let us now replace Eq. (|57p into the trajectory equa- 
tion (|56|) . Because the ground-state wave function 
does not depend on the Lagrangian trajectories, we can 
pull the functional derivative with respect to r M out of 
the quantum average. In other words, we can write 




# ) + 



SWp 



S 



V \H + W \* ) 



A. 

(58) 



where Wo is the operator 



W [g^} 



d€n(€)ff^(€,*K«K(€). (59) 



such that (tpo\W \tpo) = W . 

The functional derivative with respect to r^ at con- 
stant A. can be written as an unrestricted functional 
derivative minus counter-terms, which cancel the un- 
wanted contributions from the variation of A.. In for- 
mulas, we have 



, SA V (£) 



*o\H + W \*o) 
5 



<*o|#|*o 



(60) 



where we have used the fact that Wo does not depend on 
A.. Now we recall the identity 



j (£ ) = 



SAM') 



(61) 
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and the fact that j as well as A v [see Eq. (|4"8)) ] must 
be here evaluated at ^o- In particular, we notice that, 
according to Eq. (|48|) 



(62) 



Therefore, Eq. (|60|) can be compactly and suggestively 
restated as 



Here E c \[g^ v ] is an elastic energy defined as 



FelAZ) ■ (63) 



#el[ff^] = ("0o|#[.9a^;-4 = -Vpo]|^o) 



(64) 



where v p o = j P o/^o is the paramagnetic velocity in the 
ground-state (in a proper gauge, this quantity always 
vanishes in a one-particle system). We shall refer to F c \, 
as the elastic force acting on the systems: this is given 
by the functional derivative of the clastic energy with 
respect to the Lagrangian trajectory. 

Plugging the above definitions and formulas into the 
trajectory equation (|56[) , we finally obtain 

D t V + [(Ar) x BoL+^Vb ~ = -d^Vi . (65) 

no 

The situation is actually analogous to what happens in 
the case of a classical spring: displace it from the equi- 
librium, compute the gain in energy, get the elastic force 
from the derivative of the energy with respect to dis- 
placement. Classical elasticity theory generalizes this to 
continuum media, for which the "stress field" is the func- 
tional derivative of the energy with respect to the strain 
field.— Here, in a quantum many-body system, the cor- 
responding quantity is the functional derivative of the 
(transformed) quantum energy with respect to the tra- 
jectories of the infinitesimal volume elements. Due to 
the presence of an external static magnetic field, those 
material elements carry a non- vanishing current. 

The quantum clastic energy can be naturally split into 
a kinetic component and an electron-electron interaction 
component. In detail, we find 



(66) 



where 



1/2 



(67) 



AT?) 



2 ^J^ffi + Kfphoit, O 

{dy.\fr^){d v y/n^) , 



(68) 



7o(r, r') is the ground-state one-particle density matrix 

N 
3=2 

it is the operator of the kinematic momentum for the 
relative motion 



7T = -iV - v 0p , 



(70) 



W el [ gilv ] = \J ded€'W(r(€)-r(0)ro(€,0. 



(71) 



and 



N 

r (r,r') = N(N-l) / JJ dr 

J 3=3 



is the ground-state pair distribution function. 

The identification of the elastic energy defined by 
Eqs. (|66H72j) is the key result of this section. In com- 
parison with the non-magnetic case^^ the essential dif- 
ference is the redefinition of the kinetic contribution [see 
Eq. fl57|), Eq. flSHD, and Eq. GOD]. 



V. LINEARIZED EQUATION OF MOTION 



The equation of motion for the Lagrangian trajectory 
derived in the previous section are fully nonlinear. Here 
we present the linearized form of those equations, which 
are expected to be useful for systems performing small 
oscillations about the ground-state. 

First of all, we insert Eqs. ([46j) and (J3TJ) on the left 
hand sides of Eqs. (|55"|) and This gives 



(73) 



- n(r(t,t),t)v»(r(t,t),t)] , (74) 

where 

r(£,t)=€ + u(*,t). (75) 
Expanding to first order in u we get 

«(€,*) = «o(e)-V-[no(€)u(€)] (76) 

and 

t) = j (£) + n (£)u(£) + V x [u(£) x j (£)] . (77) 



10 



Eq. (|77[) expresses the response of the current density, 
j — jo, as the sum of two terms: the first accounts for the 
polarization and the second for the orbital magnetization 
of the quantum medium. We have thus succeeded in 
deriving the relation (J8J) between the linearized response 
of the current and the displacement field that was put 
forward in the introduction on a heuristic basis. 

We notice that, up to the linear order, there is no dif- 
ference between the Lagrangian and the Eulcrian descrip- 
tion of the displacement field, thus u(£) = u(r). In this 
spirit, we replace £ by r is all the expressions below. It 
is also worth emphasizing that Eqs. ([75)1 and (|77p are 
general and hold true independently of the elastic or any 
other approximation. 

The linearized equation of motion of the displacement 
is readily obtained 

D 2 t u{r,t) + [Au(r,t)j x B (r) + (u • V) ^ K (r) 
+ v x (u- V)Bo(r) -—F el {r,t) 



-Wi(r,t). 



(78) 



where J- e \{^, t) is the (linearized) elastic force. The terms 
on the second line result from the expansion of the clas- 
sical forces to first order in u. The linearized elastic force 
is given by 



7" e i iAl (r,t) = - J dr' 



<5 2 #el[u] 



6u,j.(r)8u v (r') 



Uv(r',t), (79) 



and can be naturally separated into kinetic and electron- 
electron interaction contributions: 



(80) 



Explicit expressions for the two components are ob- 
tained by closely following the steps outlined in Ref. 0. 
The final expressions are 

J*>) = 9 a [2Q(r) Hra (r) + Q(r)^(r)] 
- ^d v d,i [n (r)d v V ■ u(r)] 

+ l -d v {2 [V 2 n (r)] u„ M (r) + [d v n (r)} d^V ■ u(r) 

+ [0 M no(r)] 0„V • u(r) - 2d M [(d a n {r)) u„ a (r)}} (81) 

(iiv^r) = (9„u M (r) + it;, (r))/2 is the strain tensor) and 
the interaction contribution is given by 



(82) 



-^(r) = J rfr'^(r,r')K(r) - ^(r')] . 
In Eq. (jgTj) , we have defined 

^ S (r) = I(^ + ^;) 7 o(r,r')| r = r '-jV 2 no(r)V 

(83) 



3 n (r)v p o^(r)v p o^(r) , 



where 7o(r, r') is the ground-state one-particle density 
matrix [see Eq. (|69]) ]. In Eq. ([82]l 

(84) 



AV(r,r') = r (r,r')^a^(|r - r' 



where W(\r — r'|) is the interaction potential and 
r (r, r') is the ground-state pair distribution function 
[see Eq. (ff2|)]. 

Analysis of the one-particle case within the Lagrangian 
formulation is reported in the Appendix C, thus this also 
completes the analysis started in the Appendix B within 
the Euler approach. 



VI. CONCLUSIONS 

In summary, we recapitulate the essential changes that 
must be made in order to go from the quantum contin- 
uum mechanics in the absence of magnetic field to the one 
in the presence of magnetic field, within the framework 
of the "elastic approximation" : 

(i) A Lorentz-force term must be added to the equa- 
tion of motion for the displacement field, and ap- 
propriately linearized, taking into account the pres- 
ence of a non- vanishing velocity field vo(r) in the 
ground-state. 

(ii) Everywhere, the time derivative dt must be re- 
placed by the convectivc derivative D t ~ dt + 
(vo • V). The replacement must also be done within 
the Lorentz force term, compatibly with the re- 
quirements of linearization. 

(iii) The kinetic contribution to the clastic energy must 
be calculated taking into account the replacement 
of the canonical momentum operator — zV by the 
kinematic momentum tv = — iV — vo p . 

(iv) The relation between linearized current density 
and displacement field is changed from Eq. (TTJ) to 
Eq. ©. 

The Fourier transform of Eq. ([75]) yields the following 
generalized eigenvalue problem 

oj 2 u + iuj [2 (v • V) u + (u x B )] - (v ■ V) 2 u 

- [(v -Vu)xB ]-(u-V)9 M Vo(r) 

- v x (u • V) B (r) + — T cX = . (85) 

n 

Finally, we have an explicit form of all the terms: this is 
a major step forward. 

In conclusion, we have presented results that open the 
possibility to obtain the response of the current, and 
thus the excitation energies, of systems in strong mag- 
netic fields avoiding the solution of the time-dependent 
Schrodinger equation and making use only of ground- 
state properties. Given the required ground-state prop- 
erties, the complexity of the problem to be solved does 
not increase with the number of the particles in the sys- 
tem. The presented framework is expected to be useful 
in dealing with large systems and with current-carrying 
states exhibiting an clastic behavior. 
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Appendix A: Generalized eigenvalue problem 

Eq. can be represented in the form 

w 2 u-MwB-u-7C-u = 0, (Al) 
where we introduced the following notations 

u M (r) = \/no(r)u„(r) , (A2) 



Let us discuss the conditions for which the solutions 
are guaranteed to have real- valued u. For this, Eq. (jAlj) 
can be brought the quadratic form 

w 2 (u, u) + i u(u, B ■ u) - (u, K ■ u) = . (A9) 

Real- valued ui are obtained for positive definitive discrim- 
inant of Eq. (|A9[) . Since (u, B ■ u) is a purely imaginary 
quantity, the stability condition may be stated in terms 
of the stronger requirement 

(u,X:-u)>0. (A10) 

Eq. (|A10[) is not manifestly satisfied for the equations 
under considerations. Nevertheless, for Eq. ([2^]) and 
Eq. (|85|) . we expected to find stable solutions because 
those same equations are valid for small displacements 
and short-time intervals: within these conditions, the 
system must stay "close" to the the initial minimum (the 
ground-state) of the unperturbed Hamiltonian. 



B M ,,(r,r') 

*W r > r ') 



-^L=B^ v (r,r') r l_ , (A3) 



\/no(r) 
1 



\/n {r) 



l 



(A4) 



Moreover, Eq. ([85]) shares the same structure although 
(as we have explained) it is for a different dispalcement. 
The operators acting on u have properties 

S^(r,r')-6^(r,r'), B M) „(r, r') = -B^(r', r) (A5) 

and 

K^r, r') = iC*^(r, r'), fC^r, r') = t v ^{v' , r) . 

(A6) 

Eq. (|A1[) has the form of a non-standard eigenvalue prob- 
lem which may be found for example also in the analysis 
of the modes of rotating stars^ 2 - and in the analysis of the 
magneto hydrodynamic of hot plasm a 33 ' 34 . 

In terms of the scalar product [the difference with 
Eq. (j2l))) is related to the rescaling of the displacement 
field, Eq. (|A2]) ] 

(u B) ua) = J dr u* B41 (r)u A ^(r) , (A7) 
the orthogonality relation get modified as follows 

(ub, iB ■ ua) + (u>a + wb) (ub, ua) = , (A8) 

where ua and lob are assumed to be real and different 
from each. In fact, the operator defining the problem 
in Eq. (|Alj) is Hcrmitian for real ui but it is also in- 
dependent; thus, solutions of Eq. (|A2j) corresponding 
to different lu do not need to be orthogonal in the usual 
sense (as for 8 = 0). 



Appendix B: One-particle case in the Euler 
description with the standard displacement 

Here, we show that the inversion of the current re- 
sponse for one particle system can be easily worked out 
from the linearized Schroedinger equation or, equiva- 
lcntly, from the linearized Euler equations for the den- 
sities. In this way, we are able to show that the approx- 
imation put forward in the high-frequency limit provide 
the exact excitation energies for one particle systems. 

We start by observing that, the wave function may be 
written as follows 



thus, 



and 



j = riV(p + nA. 



V x v = B 



(Bl) 



(B2) 



(B3) 



As a result, the local balance equation for the linear mo- 
mentum reads as follows 

dtU + 9„ [ Vb + ~ + V ] = d t A„ (B4) 



where 



vv 



i v 2 v^ 



(B5) 



2 \M 

is the well-known Bohm potentia l 35 ' 36 . Let 

V = V , and A = A + Ai , (B6) 
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the linearization of Eq. (|B4[) yields 
dth (V-jj 



no n 



where 



v + Vl,B + v • — -«§ — 
ri n 



72 ni 



no ( ni 



(B7) 



(B8) 



(B9) 



From Eq. (|B7|) . we obtain an important relation for 
the inverse of the current-current response function: 



moreover, we may also remind that 
d t m = -V • ji . 



ji ^ 1 

no u) 
J_V 
^1 



(V-ji 



1 



2vi 



n 
Vji 



vo 
2 



■ V Jv 

2x/no 



n 



(BIO) 



This expression tells us that the frequency dependency 
includes only 1/ui and 1/tu 2 terms: therefore, it is ap- 
parent that the inversion in high-frequency limit of the 
current response function is exact for one particle sys- 
tems. 

Let us consider purely longitudinal perturbation, from 
Eq. (|B7[) we get the equation for modes described in 
terms of ji = -iunoi 



u 2 u 



V 

2 

VVi , 



(V • n u) 
1 



vo 

Yl 

2 



V (v ■ u) 

V 2 



/n Q 



- 2vt 



V • (n u) 



/no 



(Bll) 



Now, we find the exact displacement for the one-particle 
system. For this purpose, we use the expression of the 
linearized Schrodinger equation. Without loosing gener- 
ality, one may shift to zero the energy of ground-state 
energy and choose the ground-state wave function in 
Eq. (|B1[) to be a real- valued function (i.e.; (fo = 0). In 
the given gauge, the relation 



A _ Jo 

no 



vo 



(B12) 



holds true. By substitution in Eq. (jBllj) . we can verify 
that the solutions have form 



are the matrix element of the current operator evaluated 
for the eigenstates, of H and w n0 are the correspond- 
ing excitation energies. In verifying the above results, it 
may be useful to remind the continuity relations 



J On 



ZW n0 *0*n 



(B15) 



Appendix C: One-particle case in the Lagrangian 
description with the new displacement 

Let us consider the one-particle case but within the 
elastic approximation as obtained within the Lagrangian 
description. In this case, we readily arrive at 



g l/2 



^/S) ■ ( C1 ) 



It is remarkable that, Eq. (|C1[) has the same form as 
in the limit of vanishing magnetic field. Moreover, it is 
possible to directly verify that the Lagrangian with the 
elastic energy as in Eq. (jCll) : 



L = j d£n (£) 



i(Ar) 2 -(Ar)A (r)-^(r,i) 



»o 



gl/2 



no 



1/2 



(C2) 



yields the exact equation of motion for the trajectory 
r(£, t). In other words, the elastic approximation is exact 
for one- particle systems regardless the value of Bq. 

Then, the linearized equation of motion is readily ob- 
tained 

iiu + vo^dyii^ - u v d y vo^ + <9 U {vo, v u v ) 

V 2 u al3 d a u a pdi3y/no~ 



0, 



(C3) 



u 1 V0,aV0,pU a p 

u a d a V } = 9Ai !jU 



The latter expression can also be obtained from Eq. (|B4j) 
by re-expressing the response of the current in terms of 
"Lagrangian" displacement 



j\,n = noiifj, + d v (jo^Un - jo^Up) 



(C4) 



In conclusion, as expected, the same equation is obtained 
if the same displacement is employed consistently. 



J 0» 



n 



W — w n0 



where 



J On 



) 



(B13) 



(B14) 
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